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I. Introduction 

In recent years, there has been intense study of nanostructures such as 
quantum dots (QD) [1-5], where quasi- two-dimensional islands of electrons 
are laterally confined by an externally imposed potential that, in a good 
approximation, is parabolic. In Ref.[6], the electronic states of interacting 
electrons in three-electrons QD are calculated without making assumptions 
about the shape of the confining potential and dimensionality of the problem. 

Theoretical investigations of the ground states of QD have been reported 
in many papers. As for the standard Hartree and Hartree-Fock (HF) ap- 
proximations, there are doubts about their accuracy, since the exchange and 
correlation energies can be significant in QD [7,8]. 

A simple way to incorporate the interaction between electrons is to use 
the Post model [9], where inter-electron repulsion is replaced by the harmonic 
interaction [10]. For a critical analysis of this approximation, see Ref. [11]. 
The Post model [9] was used for a problem in nuclear physics in Ref. [12]. 

Numerical calculations using exact-diagonalization techniques were car- 
ried out in Refs. [13-19]. These calculations are computationally extensive 
and limited to a few (< 6) electrons. 

The ground states of an N-electron QD in magnetic fields have been 
measured up to < 50 [20]. 

The purpose of this work is to provide a rigorous lower bounds to the 
ground-state energy of N-electron QD in magnetic fields for any N. We show 
that our lower bounds for ground states agree well with the exact results of 
the diagonalization method. 

The paper is organized as follows. In Section II, we generalize a lower- 
bound method developed by Hall and Post [21] for the case of N-electron QD 
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in a magnetic field B. In Section III, lower bounds are found analytically in 
the large B limit. In Section IV, we describe our calculations. A summary 
and conclusions are given in Section V. 

II. Lower Bounds 

The Hamiltonian for iV interacting electrons confined in a parabolic QD, 
in the presence of a magnetic field B perpendicular to the dot, can be written 

as 

H = ^fl$ + \™*& E rl - y l * + E 7Y^¥T\ + 9 ^ bBS - (1) 

i=l ' 1=1 ' i<j I 1 3 I 

where m* is the electron effective mass, f2 2 = lUq + uj^/A, lu is the parabolic 
confinement frequency, uj c is the cyclotron frequency, L z is the z component 
of the total orbital momentum, e is the dielectric constant, g* is the effective 
g-factor, fit, is the Bohr magneton, and S z is the z component of the total 
spin. 

In our numerical calculations we use the effective mass m* = 0.067m e 
(m e is the free-electron mass) of GaAs QD. 

Now we introduce the center-of-mass coordinates, R — ^ J2iLi U an d 

Using 

E*? = ^ + ^E(^-^) a , (2) 

i=l iV i<j 

Etf = ^ + ^E(a-a) 2 > ( 3 ) 

1=1 j<jr 

and 

w 1 
L = J2riXpi = Rx P+— J2(ri x p< + x $ - x - x (4) 

i=l iV i<j 
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we can rewrite Eq.(l) as 

H = H cm + H re i + H z , (5) 

where the first term is the center-of-mass energy, the second term is the 
relative energy and the last term is the Zeeman energy, with H z = g*fj Jb BS z . 
H cm and H rd are given by 

P 2 m*N&R 2 m 

and 

Hrel = ^2Hij, (7) 

i<j 

where 

(Pi—Pi) 2 m*n 2 (fi — fj) 2 e 2 uj c _ _ _ _ _ _ 

(8) 

Hence we have for the ground state energy 

E = Ml + E rel + g*fi b BS z , (9) 

where 

E rel =< ij) | H rel | V >, (10) 

and ip(fi,r2, ...fjv) is the ground state wave function . Using the symmetric 
properties of ip we can rewrite Eq. (10) as 

N(N - 1) 

E rel = — <^\H 12 \ij}>. (11) 

— * 

Introducing the Jacobi coordinates Q as 

N 

3 = E^» (12) 



with 



Uij = < 



(i)- 1 if j < i + 1, 

-1 if J = i + 1, (13) 
if j > % + 1, 



we have 



2 a m*tt 2 ( 2 e 2 u C/ ? ft , 

where £i = fi — r*2. 

Projecting | ^ > on the complete basis | n >, generated by the effective 
two-body Hamiltonian H\ 2) H\i \ n >= E n \ n >, and using 

< ip | H l2 | ip >= J2 E n \< I n >< n | if; >\> E g , 

n 

we get 

E>hn+ N< " N ~ ^ E g + g*fi b BS z (15) 

where E g is the ground state energy of the effective two-body Hamiltonian 
H 12 (for the completely spin polarized states, S z = Nh/2 and E g is the energy 
of the lowest antisymmetric state of the effective two-body Hamiltonian H 12 ). 
Eq.(15) is a generalization of the Hall- Post method (which is restricted to 
the case with only interparticle forces present and no external potential) for 
obtaining lower bounds to the ground-state energy of N-electron QD in a 
magnetic field B. 

III. Large B Limit 

We introduce dimensionless units by making the following transformation: 



p = (l/a)Ci, where a = Jh/(m*uj ). 
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Using the above dimensionless notation and polar coordinates p x = p sin 6 
and p y = p cos 9, we can write the effective two-body eigenvalue problem, 

H12 | 4> g >= Eg | (f) g > as 

rr , x r 2 d2 2(£ 2 -l/4) 1 „ A 2 , 2 o7c £\, . , ~ , , 

(16) 

where A = uj c /uq, <p g (p,9) = e M u(p) / \fp, 7 C = a^Jm*c 2 / (hu )/2, and 
E = E g /(huJo). 

The two-body equation, Eq.(16), can be solved numerically to find E g for 
any arbitrary value of t. The optimal I value, restricted to odd integers for 
polarized states, minimizes the energy. The two-electron QD has been the 
subject of intensive study [3,4,8,22-26]. 

In the large magnetic field limit, i becomes large and the term 1/4 in 
(I 2 — 1/4) can be neglected [17], and Eq.(16) can be rewritten as 

{-^2+ V eff(p)Hp)=Eu(p), (17) 

where 

v <"M = w + ^ (1+ ^ )p2+2 ^- tm - (18) 

In the large- B limit the effective potential V e //, Eq.(18), has a deep minimum, 
therefore a good approximation to E can be obtained by making the Taylor 
expansion of V e ff about its minimum [17]. Thus the approximate E g is 

E g « ^[^(27eiV) 2/3 + ^IH. (19) 

Substitution Eq.(19) into Eq.(15) gives 

E > E low ^8 = ^+(iy-l)[^(2 7c iV) 2 / 3 + iv / A 2 T3]^o + ^/iB5^. (20) 
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Note that the large N limit of Ei ow is independent of magnetic field in this 
approximation (see also Ref. [16]). 

IV. Numerical Results 

We begin with the single-electron basis functions Xni(p), associated with 
Hamiltonian 

d 2 ^ 2 -l/4 1,, A 2 , 2 \i 

= ~dj? + + 4' + ~ T' (21) 

These functions were found more than seventy years ago [27], 

XM = W +1/2 e- (1/4)(1+A2/4)1/2p2 ^((l/2)(l + A 2 /4) 1/2 P 2 ), (22) 
where L l n are associated Laguerre polynomials and 

^( 1 + T>' +1 (^1 1/2 < 23 » 

In order to solve Eq.(16) we introduce 

■xSM = ^Xne((3p) (24) 
and expand u(p), Eq.(16), in the basis Xni-> l - e - we see ^ solution of the form 

M 

« M (p) = E4xf,(p). (25) 



The conventional choice for the parameter (3 is f3 = 1 (see, for example [4]). 
However, for finite M, the choice (5 — 1 is not the optimal choice. The most 
reliable (3 is obtained from 

-j- < u M I H I u M >= 0. (26) 

We apply the method, Eqs.(25) and (26), to compute the lower bounds. 



Consider a two-dimensional three-electron QD with e = 13.1 and Hujq = 
0.01 meV without a magnetic field, B=0 [18]. Let M be the number of 
functions in Eq. (25). Examples of the lower bound to the completely spin 
polarized three-electron state, S z = 3h/2, corresponding to the different M 
are given in Table I. The fast convergence is evident. Comparision of the 
converging result of Table I, E g = 0.336659 meV with exact diagonalization 
calculations of Ref.[18], E g = 0.3393 meV, shows that our lower bound is a 
very good approximation with relative error of about 0.7% for the ground 
state energy of the three-electron QD without a magnetic field. 

Now consider the GaAs QD with e = 12.4 and Tiuq = 4 meV in a strong 
magnetic field, B = 20T [19]. Examples of the lower bounds to the com- 
pletely spin polarized N-electron ground state, for up to N = 6 electrons 
are given in Table II. Numerical results Ei ow agree with large B approxima- 
tion, E (Eq.(20)) to better than 0.1%. 

From Table II, we can see that the calculated lower bounds agree well 
with exact-diagonalization results, E e ^ [19]. The relative error, A = {E e d — 
Ei ow ) / (2E e d) , is less than 2%. 

We have also calculated the chemical potential of QD, (J,a(N) = E(N + 
1) — E(N). /jLa(N) is measured by the single-electron capacitance spec- 
troscopy method [20, 28] since the transfer of the (N + l)th electron from 
the electrode to the QD occurs when the chemical potential of the electrode, 
fiE, is equal to the /ia- 

V. Summary and Conclusion 

In summary, we have generalized the Hall- Post lower-bound method [21] 
for the case of the iV-electron parabolic QD in the presence of a perpendicular 
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magnetic field. 

It is shown that our rigorous lower bounds agree well with the results 
of exact diagonalization. For example, lower bounds to the completely spin 
polarized ground state energy of the three-electron QD agree with exact- 
diagonalization results to better than 1%. For the case of six-electron QD, 
the relative error is less than 2%. 

Analytic results for the lower bounds to the ground-state energies of the 
QD in a strong magnetic field (the QD analogue of a Wigner crystal [29] 
known as electron molecule [19]) agree with numerical lower bounds to better 
than 0.1%. 
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Table I. Convergence of the method, Eqs. (25-26) for lower bounds, E g 
with increasing M for the three-electron QD with e = 13.1 and Hojq = 0.01 
meV. S z = 3h/2 and B = is assumed [18]. 



M 


1 


4 


5 


6 


7 8 


Eg, meV 


0.373368 


0.336831 


0.336681 


0.336659 


0.336659 0.336659 



Table II. Results for lower bounds Ei ow , chemical potential /ia, large B 
analytical approximation S, Eq.(20), and A = (E ed [19]—Ei ow )/(2E ed [19]) for 
N-electron GaAs QD with e = 12.4 and huj = 4 meV in a strong magnetic 
field, B = 20T. S z = Nh/2 is assumed. 



Number of electrons, iV 


E low , meV 


Ha, meV 


£, meV 


A, % 


1 


17.4810 


24.2141 


17.4810 




2 


41.6951 


28.5087 


41.6910 




3 


70.2038 


32.0032 


70.1483 


0.2 


4 


102.207 


35.1570 


102.168 


0.7 


5 


137.364 


38.2630 


137.366 


1.4 


6 


175.627 


40.6922 


175.482 


1.8 
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